% Ejemplo de solucion numerica de un sistema de ecuaciones 
% diferenciales ordinarias mediante el metodo de Runge-Kutta-Fehlberg.
% Requiere rkf45.m y arenstorf.m.

% Se ponen las condiciones iniciales
a = 0;
b = 50000;
ya = [12000;6000];
M = 10000;    %h=(b-a)/M; 1000000
tol=1e-7;
clf;
z= rk4('guanacos2',a,b,ya,M,tol); %Retorna [t, [x, y], h]
hold on;
plot(z(1, :), z(3,:),"-;Pumas;");  %y(t)
plot(z(1, :), z(2,:),"-;Guanacos;");  %x(t)
xlabel('Tiempo(dias)')
ylabel('Unidades de batalla')
title('Batalla de Iwo Jima')
hold off;

